Exotic fission properties of highly neutron-rich Uranium isotopes 
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The series of Uranium isotopes with N = 154 ~ 172 around the magic number N=162/164 are 
identified to be thermally fissile. The thermal neutron fission of a typical representative 249 U of this 
region amenable to synthesis in the radioactive ion beam facilities is considered here. Semiempirical 
study of fission barrier height and width shows this nucleus to be infinitely stable against spontaneous 
fission due to increase in barrier width arising out of excess neutrons. Calculation of probability of 
fragment mass yields and microscopic study in relativistic mean field theory, show this nucleus to 
undergo a new mode of thermal fission decay termed multifragmentation fission where a number of 
prompt scission neutrons are simultaneously released along with the two heavy fission fragments. 

PACS numbers: 21.10.Dr, 24.75. +i, 25.85.W, 21.60.-n, 24.10.Cn 



Exploration of nuclear landscape with radioactive ion 
beam facilities is currently underway in several laborato- 
ries around the world. Presently about 2000 nuclei are 
known, and different mass models predict the survival of 
another 5000 nuclei with varying life time which could 
be synthesised in the laboratory [l[. Most of the nuclei 
will occur in the neutron-rich side of the nuclear chart. 
Many of these nuclei will have unusual neutron to proton 
ratio, quite different from the usual ones in the valley 
of stability. Therefore with the synthesis of such nuclei 
one may expect to observe new nuclear phenomena and 
novel features of nuclear dynamics. The two Uranium 
isotopes 233 U and 235 U and Plutonium isotope 239 Pu in 
the actinide region are suitable for energy production be- 
ing thermally fissile in nature. Here we investigate to see 
if heavier neutron-rich isotopes of Uranium could exist 
having thermally fissile properties, and if so what would 
be their fission decay properties. 

On the basis of fission barrier Bf and neutron sep- 
aration energy S n systematics, we find that the chain 
of Uranium isotopes with neutron number N=154 ^172 
possess the thermally fissile property 0]. These isotopes 
span on either side of N=162/164 which has been pre- 
dicted to be magic in numerous theoretical studies car- 
ried over the years [l| . The effect of close shell manifests 
in relatively higher values of S n and lower values of Bf 
(favouring more bound compact configuration) rendering 
the system thermally fissile, which may be the case with 
the above series of Uranium isotopes. All these nuclei 
are stable against alpha decay and several of them in the 
lower side have beta-decay half-life of several tens of sec- 
onds [4J. We have chosen 249 U, which leads to 250 U with 
capture of thermal neutron for the present study. This 
isotope is only 11 neutrons away from the naturally oc- 
curing known Uranium isotope 238 U, and therefore likely 
to be produced with RIB facilities under construction. 
The question of how the neutron-rich isotopes of the ac- 
tinide nuclei will decay by fission has not been addressed 
before. The neutron separation energy S n provides a 
measure of excitation of the compound nucleus. In our 
study we use the fission barrier from the extensive study 
of Howard and Moller (HM) [2| which is being widely 



TABLE I: Fission barrier Bf (in MeV) obtained in finite- 
range droplet model (FRDM) [5J and Howard and Moller 
(H.M.) 0] are compared with experimental data [fj. 



nuclei 


B f (expt.) B f 


(H.M.) Bf 


(FRDM) 


228^ 


6.50 




7.43 


230^ 


7.00 


7.79 


7.57 


232^ 


6.30 


7.40 


7.63 


234 Th 


6.65 


6.83 


7.44 


232 yj 


5.40 


6.31 


6.61 


234TJ 


5.80 


6.03 


6.79 


236yj 


5.75 


5.74 


6.65 


238 yj 


5.90 


5.83 


4.89 


240tj 


5.80 


5.92 


5.59 


238 p u 


5.30 


5.25 


4.85 


240 Pu 


5.50 


5.92 


4.74 


242 p u 


5.50 


5.48 


5.25 


244 Pu 


5.30 


5.31 


5.78 


246 p u 


5.30 


5.06 


6.27 


242 Cm 


5.00 


5.56 


4.24 


244 Cm 


5.00 


5.56 


5.05 


246 Cm 


4.70 


5.40 


5.69 


248 Cm 


5.00 


5.08 


6.07 


250 Cm 


4.40 


4.53 


5.51 



used in literature. It is also the only exhaustive calcula- 
tion which gives barrier for the whole range of known and 
unknown nuclei where fission decay can occur. We have 
compared in Table I the HM values of Bf with experi- 
ment for 18 actinide nuclei of the four elements Th, U, Pu 
and Cm for which data are available Q. These data are 
for the external barriers which in general correspond to 
the maximum barrier height calculated by Howard and 
Moller. We have also included in the table the Bf val- 
ues calculated by Moller et al [5j in finite-range droplet 
model (FRDM) for a relative comparison. It is interest- 
ing to see that HM values are in much better agreement 
with experiment than those of FRDM, although the lat- 
ter are comparatively more recent. For the Uranium iso- 
topes which is of special interest in the present study, the 
agreement is remarkably satisfactory. Therefore we have 
used them in Fig. 1, shown as closed circles for all the 
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Uranium isotopes with the neutron number in the range 
140 ~ 180. For the neutron separation energies SL, we 
have used the three mass formulae HM Q , FRDM [3] and 
infinite nuclear matter (INM) || . It is well known that if 
S n < Bf, then the nucleus cannot undergo thermal neu- 
tron fission. The fission threshold E nm = Bf — S n has 
to be overcome by impinging with an energetic neutron, 
in order that the nucleus will undergo fission. However, 
if S n > Bf, then thermal neutron (with practically zero 
energy) can cause fission. It is interesting to see in Fig. 1 
that, the isotopes of Uranium in the range N — 154—172 
are thermally fissile, a feature emanating from the close 
shell nature of N= 162/ 164. To show that this feature is 
not accidental or specific to Uranium, but a general one 
arising out of the close shell structure of N=162 (or 164), 
we have presented the case of Th-isotopes in the same 
figure. 

We now consider the nature of the fission decay mode 
of 250 U, which is primarily governed by the profile of 
the fission barrier. The height and width of the fis- 
sion barrier which is supposed to be parabolic in nature 
have to be obtained Q- We have followed a semiempir- 
ical method to get the width of the barrier from the 
systematics of the known experimental fission half-lives, 
and extrapolated them to neutron-rich region of inter- 
est. In spontaneous fission the fissioning fragment will 
see the maximum barrier which will determine the tun- 
neling probability more decisively. Hence, the maximum 
barrier given byHM can be considered as the effective 
parabolic barrier governing the decay, which should be 
appropriate for the present study. The fission half-lives 
can be calculated using the relation [?J tx/2 = ln2/np, 
where n is the number of barrier assaults by the decay- 
ing fragment related to the barrier curvature energy Huj 
by nh = Huj/2tt, and p is the penetrability of the barrier 
given by p = [1 + exp{2irB j / 'Hui}^ 1 . It may be noted 
that Huj is a measure of width of the barrier; smaller 
the Huj larger is the barrier width, and hereafter we refer 
it as width to bring out a physical picture. Taking the 
HM values [|J of Bf and using the experimental Tx/% @, 
we get the values of Huj. The systematics of Huj versus 
Bf so obtained are shown in Figure 2 for various even- 
even actinide nuclei. It is indeed quite revealing that the 
plot shows a linear behaviour of Huj with Bf for a given 
Z, with progressively increasing slope with the increase 
of proton number of the elements from 90 to 96. For 
the next element Cf with Z=98, the linear behavior gets 
fuzzy. However, the mean follows the trend with a higher 
inclination. The trend is more conspicuously restored for 
Fm with Z=100. This deviation coincidentally correlates 
well with the fission mass yield systematics, where con- 
siderable deviation from the standard well-defined two 
peaks occurs for Fm isotopes. The width Huj for any 
isotope with calculated fission barrier may be obtained 
by extrapolation of the linear graph. Since 250 U has 158 
neutrons and of considerable distance from the close shell 
N=162 (or 164) and itself only 12 neutrons away from the 
known 238 U this extrapolation should be quite meaning- 



ful and reliable. For 250 U with a fission barrier of 4.3 
MeV |2[, we obtained the value of Huj as 0.225 MeV from 
Fig. 2. Now we can construct the parabolic barrier of 
base width Ar somewhat schematically, using the value 
of fojj following the relation 1/hco = d 2 V/dr 2 , where V 
is the potential energy. The barrier so obtained is rel- 
atively flat and wide compared to the fission barrier of 
236 U with Huj = 0.357 MeV, obtained from the experi- 
mental half- life Q and Bf = 5.74 MeV from Howard and 
Moller 0. The normal nuclei like 2 30,232 Th or 235,238tj 
have Huj k. 0.4 — 0.5McV. Thus, the neutron-rich heavy 
isotope 250 U has considerably lower value of Huj, and con- 
sequently relatively larger width, arising out of excess of 
neutrons. This flattening of fission barrier makes the nu- 
cleus stable against spontaneous fission decay, because of 
decreasing penetrability. The spontaneous fission decay 
half- life calculated with these height and width comes out 
to be 5.7 x 10 24 years, which is several order of magnitude 
higher than that of 236 U. However, due to the decreased 
fission barrier of 4.3 MeV for 250 U, with a minute induce- 
ment by a thermal neutron, fission decay will occur. This 
exotic feature is due to excess of neutrons, which is not 
the case in the normal nuclei in the valley of stability. 

Fission studies 0, [l(| of 235 U with thermal neutrons 
have shown that neutron emission from the neck region 
is many times (order of magnitude) larger than the alpha 
particle emission. Since clustering probability increases 
at low density, alpha particle emission itself is larger com- 
pared to that of proton. Therefore, it is believed that 
neutron-rich neck is produced during the scission. This 
picture is further reinforced by the polarization effect in- 
duced by the Coulomb repulsion between the two newly 
formed fragment nuclei. It is also supported in micro- 
scopic studies throu gh o ur calculation in relativistic mean 
field (RMF) theory [l2| as will be shown aposteriori. 
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FIG. 1: Fission barrier Bf, and binding energy of the last 
neutron Si n as a function of mass number A for Th and U 
isotopes. The Bj are taken from || and Sin are taken from 
Refs. 0, [|, and [H for HM, INM, FRDM and RMF 
model respectively. 



3 




FIG. 2: The barrier curvature energy hui as function of fission 
barrier Bf. The Bf values and the experimental half-lives 
taken from Howard and Moller and Ref. [IJ, respectively. 
The star on the line for Uranium in figure, denotes the value 
of Huj to be 0.225 MeV for its B f = 4.3 MeV corresponding 
to 250 U. 
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FIG. 3: Q/ iS3 . -value distribution for Jb U and U as a func- 
tion of A\ fragment in the binary decay A — > Ai + A2. The 
binding energy used for calculation of Q— value is taken from 
[U, H and Refs. [| for RMF, FRDM and INM, respec- 
tively. The vertical line marks the neutron drip-line for the 
corresponding element in each panel. 



As is well known, the major driving force for the de- 
cay of a nucleus is the Q— value of the reaction. The 
probability of fragment mass yield in a given channel 
is directly related to the Q— value. We, therefore, cal- 
culate here the Q— value systematics of the fission of a 
nucleus (A,Z) decaying to (A\,Z\) and (A 2 ,Z 2 ) defined 
as Qf(A, Z) = BE{A U Z X ) + BE(A 2 , Z 2 ) - BE(A, Z). 
In Figure 3, we have plotted the Q— values of the binary 
decay into two fragments A\ and A%, as a function of 
the mass number of the A\ fragment, for all the relevant 
elements with even values of Z\, starting from 34 to 46. 



FIG. 4: The matter density distributions p = p„ + p p (in 
fm -3 ) and a = p n /p P of 250 U obtained in the RMF calcula- 
tions with NL3 force. 



The complimentary fragment (Z 2 , A 2 ) is thereby fixed. 
Since the yield falls rapidly with the decrease in Q— value 
for an element, we have only shown the distribution of 
Q— values lying above 90% of the highest values. For 
the sake of comparison, the Q— value distributions for 
both 236 U and 250 U are shown in Figure 3. To ensure 
that our conclusions remain quite general and valid, and 
independent of any specific mass formula, we have used 
the masses predicted in the three mass models RMF [ll[ , 
FRDM and INM. The Q-value distributions for 236,250u 
are presented in Figure 3 as dashed, solid and dotted 
curves for RMF, FRDM and INM mass formulae respec- 
tively. In the figure, three corresponding vertical lines 
mark the drip-lines of the respective elements predicted 
in the three mass models. It can be seen that the drip 
lines for all the three mass models agree within ~ 3 — 5 
neutrons, except in case of Z=40 isotopic chain for the 
FRDM mass model. In most of the cases, the drip lines 
fall inside the Q— value distributions and in some cases 
they touch the outer fringe shown in the figure with the 
exception of Z=34 and 40 where the FRDM drip lines 
are somewhat away. Thus all the isotopes lying to the 
right of the drip lines will be unstable against sponta- 
neous release of neutrons from the fragments at scission 
which is not the case with 236 Fj for which the drip lines 
are far away (See Fig. 3). In the usual fission process 
of 236 U, neutrons are emitted from the fragments after 
they are fully accelerated. But in the present case of 
250 U, a certain number of neutrons will be simultane- 
ously produced along with the two heavy fragments sig- 
naling a new mode of fission decay which may be termed 
as multi-fragmentation fission. An order of magnitude of 
these prompt multi-fragmentation neutrons as estimated 
from the mass yield plot for 250 U shown in Figure 3, turns 
out to be about 2 to 3 neutrons per fission. These are 
the additional neutrons apart from the normal multiplic- 
ity of neutrons emitted from the fragments. In case of 
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neutron-induced fission of normal 233 U and 235 U nuclei, 
the fission neutron multiplicities are of the order 2.5 @. 
This number is, therefore, more than doubled in case of 
250 U fission, which will have important implications on 
the energetics of the fission process. This phenomenon 
will be more prominent in heavier Uranium isotopes due 
to availability of more neutrons. Although, the beta de- 
cay life-time of 250 U is few tens of seconds [4( it is much 
larger than the nucleon decay life time, which is of the 
order of 10 -17 seconds and therefore it will have an im- 
plication in the r— process nuclear synthesis and conse- 
quently stellar evolution. 

We used RMF theory with NL3 interaction to study 
the evoluation of density as the nucleus undergoes dis- 
torsion in its shape on its path to fission, to see if micro- 
scopic calculation supports the above picture. Recently, 
calculations on nuclear densities obtained in the RMF 
studies have yielded useful insight and results on nuclear 
structure and dynamics [12|, which we follow here. With 
this view we carried out RMF calculations for succes- 
sively increasing deformation /3 2 starting from the ground 
state. Now we present our results of such calculations on 
the total (neutrons + protons) matter density distribu- 
tions p = p n + p p , and the ratio of neutron to proton 
density variable a — p n /pp as function of the deforma- 
tion parameter /3 2 . In Figures 4, we have presented the 
results of such calculations for 250 U, with specified col- 
umn on the right of the distributions showing the scale 
through different colours. From the figure it is clear 
that 250 U gets more and more elongated as in the usual 
liquid drop picture of fission and finally splits into two 
parts with a neck connecting them for the deformation 
/?2 = 7.4. It is interesting to examine the composition 
of the neck in this configuration. The density of the 
neck obtained in our calculation and also evident from 
the picture is p = 0.02383/to~ 3 which is quite low as 
expected. And the neutron to proton density ratio is ob- 
tained as Pn/Pp — 3.54703 which can be contrasted with 
the average density ratio of ~ 1.54 in its ground state. 
Thus these microscopic studies coroborate the neutron 



rich-ness of neck as found in experiment and other stud- 
ies [H, [l(|. Similar studies for 236 U gives the density of 
the neck p = 0.02440/to -3 and the neutron-proton den- 
sity ratio p n /pp — 2.73089. In case of 250 U, the neck is 
much richer in neutron since 14 extra neutrons are avail- 
able compared to 236 U which is in conformity with expec- 
tation. This favours simultaneous emission of neutrons 
along with the two fragments at scission, strongly sup- 
porting the multifragmentation fission process predicted 
above. 

In conclusion, we have identified the chain of Uranium 
isotopes with neutron numbers N=154 to 172 which are 
thermally fissile. This is a likely manifestation of the 
close shell nature of the magic number 162/164. We 
have chosen 249 U as a representative nucleus to study 
the thermal fission decay. The fission decay properties 
of such neutron-rich nuclei, in particular in the actinidc 
region away from the valley of stability, have not been ad- 
dressed before, which has become important in the con- 
text of RIB programs. Its fission barrier profile has been 
shown to be relatively flat and wide compared to 236 U 
yielding a half-life of 5.7 x 10 24 years, which makes it 
extremely stable against spontaneous fission, while ren- 
dering it highly vulnerable to thermal neutron fission— a 
unique property indeed. On the basis of the probability of 
mass yield and microscopic RMF calculations, strong ev- 
idence of a new mode of fission decay is revealed, where in 
addition to two heavy fragments, 2 to 3 scission neutrons 
will be simultaneously emitted, which may be termed as 
multifragmentation fission. These extra neutrons are in 
addition to the normal multiplicity of neutrons emitted 
by the excited fission fragments which is of the order of 
2.3 to 2.5. Thus the doubling of the neutron emission 
per fission will have implications for the r— process nu- 
cleosynthesis in the steller evolution. Whether such a 
nucleus presents an attractive possibility as a source of 
energy production needs to be examined. The above phe- 
nomenon is a general one, not restricted to 250 U and is 
likely to be more prominent in heavier Uranium isotopes. 
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